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Abstract 

We describe the double wall tube with cylindrical dislocation in the framework 
of the geometric theory of defects. The induced metric is found. The dispersion 
relation is obtained for the propagation of torsional elastic waves in the double wall 
tube. 


1 Introduction 

Ideal crystals are absent in nature, and most of their physical properties, such as plasticity, 
melting, growth, etc., are defined by defects of the crystalline structure. Therefore, a study 
of defects is a topical scientific question of importance for applications in the first place. 
At present, a fundamental theory of defects is absent in spite of the existence of dozens 
of monographs and thousands of articles. 

One of the most promising approaches to the theory of defects is based on Riemann- 
Cartan geometry, which involves nontrivial metric and torsion. In this approach, a crystal 
is considered as a continuous elastic medium with a spin structure. If the displacement 
vector field is a smooth function, then there are only elastic stresses corresponding to 
diffeomorphisms of the Euclidean space. If the displacement vector held has disconti¬ 
nuities, then we are saying that there are defects in the clastic structure. Defects in the 
clastic structure are called dislocations and lead to the appearance of nontrivial geometry. 
Precisely, they correspond to a nonzero torsion tensor, equal to the surface density of the 
Burgers vector. Defects in the spin structure are called disclinations. They correspond to 
nonzero curvature tensor, curvature tensor being the surface density of the Frank vector. 

The idea to relate torsion to dislocations appeared in the 1950s [1-4], This approach 
is still being successfully developed (note reviews [5-11]), and is often called the gauge 
theory of dislocations. 

Some time ago we proposed the geometrical theory of defects [12-14] . Our approach is 
essentially different from others in two respects. Firstly, we do not have the displacement 
and rotational angle vector fields as independent variables because, in general, they are not 
continuous. Instead, the triad field and SO(3)-connection are considered as independent 
variables. If defects are absent, then the triad and SO(3)-connection reduce to partial 
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derivatives of the displacement and rotational angle vector fields. In this case the latter 
can be reconstructed. Secondly, the set of equilibrium equations is different. We proposed 
purely geometric set which coincides with that of Euclidean three dimensional gravity with 
torsion. The nonlinear elasticity equations and principal chiral §0(3) model for the spin 
structure enter the model through the elastic and Lorentz gauge conditions [14-16] which 
allow to reconstruct the displacement and rotational angle vector fields in the absence of 
dislocations in full agreement with classical models. 

The advantage of the geometric theory of defects is that it allows one to describe single 
defects as well as their continuous distributions. 

In the present paper, we consider propagation of torsional elastic waves in double wall 
tube with cylindrical dislocation. This defect was first described in P2|. The Schrodinger 
equation for the double wall tube was solved in [18] and applied to double wall nanotubes. 
A similar problem was also solved for the cylindrical waveguide with wedge dislocation 

ESI- 

1.1 Double wall tube 

Let us describe double wall tube with cylindrical dislocation in the framework of the 
geometric theory of defects. 

We consider cylindrical coordinates = {r,ip,z}, /i = 1,2,3 in tree dimensional 
Euclidean space M 3 . Let there be two thick tubes ro < r < rq and r 2 < r < rq of elastic 
media, each axis coinciding with the z axis. We suppose that r 0 < rq < r 2 < r 3 (see 
Fig(T|a, where a section z = const is shown). Now we make one tube with the inside 




Figure 1: Section z = const of double wall tube before (a) and after (b) defect creation. 

cylindrical dislocation in the following manner. We stretch symmetrically the inner tube 
and compress the outer one. Then glue together the external surface of the inner tube with 
the internal surface of the outer tube. Afterwards the media comes to some equilibrium 
state. Due to rotational and translational symmetry we obtain one tube r in < r < r ex 
with the axis which coincides with the z axis (see FigjTj, 6). Radii of cylinders constituting 
tube surfaces are mapped as follows 

A) •-> An, r i, r 2 | —t A, r 3 Ax- 

The gluing is performed along the cylinder r*, and there is cylindrical defect (dislocation) 
because part of the media between tubes is removed. 
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The obtained double wall tube with cylindrical dislocation is rotationally and trans- 
lationally symmetric. 

The constructed model of the tube with cylindrical dislocation can be considered as 
continuous model of double wall nanotube (for a general review, see (2DJ12T, 23j). Consider 
double wall nanotube having two atomic layers. Suppose the inner layer has 18 and outer 
layer has 20 atoms which are shown in FigjTJfr by points. Natural length measure here is 
the interatomic distance. Then the length of a circle has a jump when one goes from inner 
to outer layer. In the geometric theory of defects, it means that the metric component g w 
is not continuous in cylindrical coordinates. The corresponding model will be described 
below. 

To find radii r in , r*, and r ex we have to solve the classical elasticity problem. 

Let us define the displacement vector field by u l (x), i — 1, 2, 3, 

y l i-> x l — y l + u l (x), (1) 

where y l and x l are coordinates of a point before and after deformation respectively. We 
consider the displacement field as a vector function on points of media after deformation 
and gluing. This is more adequate because the resulting media after gluing is a connected 
manifold (before the gluing procedure, each tube represents a connected component). In 
equilibrium state, the vector displacement field satisfies equation 

(1 - 2a) Am + didjV? = 0, (2) 

where a is the Poisson ratio and A is the Laplacian. For convenience, we consider com¬ 
ponents of the displacement vector field with respect to the orthonormal basis 

u = u r e? + + u z e z , 


where 


Cr &(n fA, 


We denote indices with respect to the orthonormal basis by hat: 


= M = {r,<p,z}. 


The Latin indices referred to an orthonormal basis are rased and lowered by Kronecker 
symbol: m := vJSji. 

The divergence and Laplacian have the following form in cylindrical coordinates 
ViU 1 = - d r {ru r ) + -d v u^ + d z u z : 


A ztf- d r (rd r Ur) T T d zz Ur 2^^ 

1 1 12 

A u z = ^d r (rd r u z ) + ^d^u z + d 2 zz u z . 


(3) 


From the symmetry of the problem, we deduce that only radial component of the dis¬ 
placement field differs from zero, and it does not depend on the angle ip and z coordinates: 


{u 1 } = {u r := u(r),u* = 0 ,u z = 0}. 
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Equation (|2]) for zero u$ and uz components are automatically satisfied. It is easy to 
check that the radial derivative of the divergence, 


1„, s\ „ 2 1„ 1 


drdjU J = d r -d r (ru ) = df. r u + -d r u -- -u, 


coincides with the Laplacian 


A Uf = -d r (rd r u) - -u = dl u H— d r u - -u. 


Therefore the radial component of Eq.([2]) takes the form 


d r ( - d r (ru ) ) = 0 . 


(4) 


A general solution of this equation depends on two integration constants: 

C2 

u = C\r H-, c 19 = const. 

r 

Note that the equilibrium equation Q does not depend on the Poisson ratio a. This 
means that the cylindrical dislocation is the geometrical defect. 

Boundary conditions have to be imposed to fix the integration constants. Let us 
introduce notation for inner and outer tubes: 


u = 



r-m <r < r*, 
r* < r < r ex . 


Now boundary conditions are to be imposed. We assume that the surface of two wall 
nanotube is free, i.e. the deformation tensor is zero on the boundary: 


du;. 


dr 


= 0 , 


du P 


dr 


= 0 . 


(5) 


We assume also that the media is in equilibrium. It means that elastic forces on the gluing 
surface must be equal to zero: 


du\ r 


dr 


du P 


dr 


( 6 ) 


Each of boundary conditions ([5]) defines one integration constant for internal and 
external solutions: 


Min 


a 



> 0 , 


Mex — 



< 0 , 


a = const > 0, 
b = const > 0. 


(7) 


Signs of the integration constants a and b are chosen in such a way that displacement 
vector is positive and negative for inner and outer tubes respectively. This in agreement 
with the imposed problem. 
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Substitution of obtained solutions (JT]) into the gluing condition (J6]) defines the ratio of 
integration constants: 


A = A 


arf n + b 
' ar 2 + b 




b = ar, 


r 2 — r 2 


ex ,v.2 


The entirety condition for the media is 


A = a + a ( r* + ^ ] , 


, , 1 a 

A = r 2 - b[ - b — 


( 8 ) 


( 9 ) 


These equations allow to find the distance between initial tubes which characterize the 
cylindrical dislocation: 

l:=r 2 -r 1 = 2ar /“~_ r % (10) 

Ax A 

where expression for b ([8]) is used. Afterwards we find the integration constants: 


a = 


l 

2 r* 


r 2 _ 2 
'ex ' * 

9 2 ’ 

(pi /p" 

'ex 'in 


b 


/ t *2 - r *2 — ^,2 

<,/ ex ' * 'in 

2 r* r 2 x — rf n ' 


Thus we find the displacement vector held 


u(r) 


a (r + ^) >0, 

-"(8 + i ) <0 - 


An < r < r*, 
r* < r < r ex 


( 11 ) 


( 12 ) 


for double wall tube where constants a and b are given by Eqs. 0- Qualitative behaviour 
of this vector held is shown in Fig|2ja. Differentiation of this vector held in domains 
ri n < r < r*, r* < r < r ex and its extension to the point r* by continuity yields the 
function 


du _ \ a ( X ~ > °’ An < r < r m , 

dr A < r < r ex , 


(13) 


which is depicted in Figj2|&. 

Note that media entirety condition 0 leads to the jump of the displacement vector 
held at the point r*: 

l ■= r 2 - n = u in (n) - u ex (n). 


Since vector held u has discontinuity at the point r*, the formal derivative of u contains 
S(r — a). This (5-function is thrown away in the geometric theory of defects [Q]. 

Double wall tube with cylindrical dislocation is parameterized by four constants r 0 , 
r h r 2 , A or r in , r*, r ex , l. Formulae 0 dehne constants a, b and, consequently, the 
displacement vector held through the second set of parameters. It follows from definition 
(|TJ) that there is one-to-one correspondence between two sets of parameters. 

Now we calculate the metric induced in double wall tube. By definition, it has the 
form 

( h^ x ) = |( 14 ) 
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Figure 2: Qualitative behaviour of the radial component of the displacement vector held 
for the double wall tube (a). The derivative of the displacement vector held (6). 


where g pa (y ) is the Euclidean metric in cylindrical coordinates. The relation between 
coordinates before and after defect creation is given by equality (Jl|. Using explicit form 
of displacement vector held (12) we find the induced metric in the double wall tube 

ds 2 = (1 — v) 2 dr 2 + (r — u) 2 d(p 2 + dz 2 . (15) 

= (1 — v) 2 of this metric is continuous function but its derivative has 


The component g 
a jump at r*. The component g w = 


u 


is discontinuous at the point r = r*. 


The volum element for metric (15) is 


VW\ = (1 - v)(r - u ). 


The right hand side of this relation is positive because both multipliers are positive. The 
second multiplier r — u = y > Ois positive by construction. The hrst multiplier is also 
positive. Indeed, the function v has the maximum at r = r*. At this point, the following 
inequality holds 




< 1 , 


where expression (J9| for a is used. 

The circumference is a geometric invariant. It is equal to 27r(r — w(r)) for metric (14). 
When we go from the inner tube to the outer one it has the jump 2nl where l is the 
distance between tubes before the dislocation is made. This observation agrees with the 
continuous model of double wall tube. 


1.2 Torsional waves 

Here we consider torsional waves in the double wall tube with cylindrical dislocation 
described in the previous section. We denote the displacement vector held by the new 
letter w because the letter u was used in the previous section for the displacements 
corresponding to the defect creation. The total displacement vector held is equal to the 
sum u + w where u corresponds to the defect creation and w describes oscillations in the 
double wall tube with cylindrical dislocation. By definition, the displacement vector held 
w satishes the wave equation 

p 0 Wi - pAwi - (A + p)ViVjW j = 0, (16) 
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where the covariant derivative V* := e^j V M and Laplace-Beltrami operator A := g fJ,u 'V l _ l 'V u 


are defined by metric (15) of double wall tube. 


This equation is covariant with respect to changing of coordinate systems. It can be 


solved in cylindrical coordinates r, (p, z after the defect creation with metric (15). However 


it is easier to follow another way. We solve the wave equation in cylindrical coordinates 
y, ip, z where y denotes the old radial coordinate before the defect creation 

y := r + u, 

and afterwards impose necessary boundary conditions. It is easier because the metric is 
Euclidean in the initial coordinate system. 

Let us consider torsional waves. In this case only angular component of the displace¬ 
ment vector field differs from zero: 

{w 1 } = {w r = 0 , wA = w^(t,r),w z = 0 }. 

From symmetry consideration, the angular component w^ does not depend on p and z. 


For this vector field the f and z components of Eq.(16) are automatically satisfied. It 


follows that the dilation for torsional waves is equal to zero 

e := diW l = 0 , 

and consequently torsional waves take place without media compression. 


We look for solution of Eq.(16) in the plain wave form 


wp= re [W(y)e< kz ~^], (17) 

where W (■ y) is the amplitude, k G M is the wave vector, and w € Mis the frequency of the 


wave. Then wave equation (16) in cylindrical coordinates reduces to the Bessel equation 


d 2 U dU , 9 0 

(18) 

dr * +T dr + (K r l)U ~ °’ 

2 7 2 2 d 1 

K := —— k , c T := —. 

c t po 

(19) 


where 


A general solution of this equation depends on two integration constants. Therefore 
general solutions for inner and outer tubes are 


W = 


Win = CMw) + C 2 N l (ny), r 0 <y < n, 
W ex = C 3 Ji(ny) + CUNi(k,ij), r 2 < y < r 3 , 


( 20 ) 


where J\ is the Bessel function of the first kind and first order, N\ is the Neumann function 
of the first order (see, i.e. |23]), and 61 , 2 , 3,4 are integration constants. 

To find the integration constants we impose boundary conditions. The boundary 
surfaces are assumed to be free, i.e. the deformation tensor on the boundary must be zero 


dW w 


dr 


= 0, 


dW P . 


dr 


= 0 . 


Because 


dW _ dy dW 
dr dr dy 


= (1 


dW 
' dy ’ 
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and r;(r; n ) = u(r ex ) = 0, these equalities in the initial coordinates take the form 


CMzo) + C 2 N[(z 0 ) = 0 
C 3 J[(z 3 ) + C±N[(zz) = 0, 


where 


* := Ky, 

and prime denotes differentiation with respect to the argument z. These equalities define 
two integration constants: 


C'2 — —koCi, fco 

C 4 = ~k 3 Ci, k 3 


J'l(z 0) 
^(*b)’ 
J[(z 3 ) 
N[(z 3 y 


( 21 ) 


On the gluing surface, we impose two conditions: entirety and equality of stresses, 


W in (n) = W ex (r *), 


dW u 


dr 


dW P ~ 


dr 


As a result we get two equations 


[Ji(^i) — koNi(zi)] — C 3 \J\{z-i) — k 3 Ni(z 2 )~\ — 0, 

cip^ij-fcoJvjM] =0. 


( 22 ) 


(23) 


The necessary and sufficient condition for this system to have a nontrivial solution is the 
equality of its determinant to zero: 


[Ji{z 1 ) ~ k 0 Ni(zi)] [,J[(z 2 ) - k 3 N[(z 2 )\- 

- [Ji(z 2 ) - k 3 Ni(z 2 )] [J[(zi) - koN[( Zl )] = 0. (24) 


For given parameters of the double wall tube r 0 , rq, r 2 r 3 , the obtained relation is the 


equation for the constant n. Let n be a root of Eq. (24), then the equality 


u 


= c T \/k 2 + K 2 . 


(25) 


dehnes the dispersion relation. 

The phase velocity of torsional waves v := ui/k is easily found from dispersion relation 

(|25b: 


v = c T 

The group velocity is also easily found 


i + £. 


du 

s dk v 


(26) 


(27) 


We see that the phase velocity of torsional waves is always greater then the velocity of 
transverse waves, and group velocity is smaller. Dispersion relation (25) does depend on 
double wall parameters through Eq.(24). 














2 Conclusion 


We found the induced metric in double wall tube with cylindrical dislocation in the 
framework of the geometrical theory of defect. Though components of this metric are 
not continuous functions, the three dimensional Einstein equations are well defined am 
Afterwards propagation of torsional waves in double wall tube is described. The presence 
of the cylindrical dislocation inside the double wall tube leads to changing of the dispersion 
relation. 

Double wall tube may be useful as a continuous model of double wall nanotubes. 

This work was supported by the Russian Science Foundation (project 14-11-00687) in 
Steklov Mathematical Institute. 
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